Effect of the quantistic zero— point atomic motion on the opto— electronic properties of 

diamond and trans— polyacetylene 
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The quantistic zero-point motion of the carbon atoms is shown to induce strong effects on the 
opto-electronic properties of diamond and irans-polyacetylene, a conjugated polymer. By using an 
ab initio approach, we interpret the sub-gap states experimentally observed in diamond in terms 
of entangled electron-phonon states. These states also appear in trans-poly acetylene causing the 
formation of strong structures in the band-structure that even call into question the accuracy of 
the band theory. This imposes a critical revision of the results obtained for carbon-based nano- 
structures by assuming the atoms frozen in their equilibrium positions. 
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Carbon-based nano-structures represent the natural 
candidates to replace silicon-based materials in the de- 
vise of efficient opto-electronic devices. Conjugated poly- 
mers, in particular, have been shown to have peculiar 
properties related to the very fast relaxation of charge 
carriers due to the electron-phonon coupling 1]. This 
ultra-fast dynamics has made possible to use conjugated 
polymers to devise efficient light emitting diodes [2| or 
nano-scale optical switches [3| . 

Despite the rapid development of technological ap- 
plications the role of atomic vibrations in carbon- 
based nano-structures has been only treated in a semi- 
empirical manner, boosted by the essential question of 
the mobility of charged carriers in organic devices 
However, these approaches are based on Hamiltonians 
that rely on parameters which are difficult to extract from 
experiments and clear-cut conclusions are still elusive. 
In contrast, the most accurate, parameter-free and up- 
to-date description of the electronic properties of bulk 
and nano-sized materials, is based on ab initio methods. 
These techniques benefit of the predictivity and accuracy 
of Density-Functional Theory (DFT) [5j merged with 
Many-Body Perturbation Theory (MBPT) @. The goal 
of the ab initio methods is to describe and predict in 
a quantitative manner, the opto-electronic properties of 
any electronic system, starting from its atomic configu- 
ration. The result is a wealth of techniques like the GW 
method H that has been successfully applied to a large 
number of different systems, among which carbon-based 
nano-structures. 

In the GW approach, as well as in other applications 
of the ab initio methods, a standard approximation is 
to assume the atoms frozen in their equilibrium posi- 
tions. Many years ago Q, however, the pioneering works 
of Heine, Allen and Cardona (HAC) pointed to the fact 
that, even when the temperature vanishes, the quantis- 
tic zero-point motion of the atoms (the zero-point mo- 



tion effect) can induce large corrections to the electronic 
levels, making purely electronic theories (like the GW 
method) inadequate. Nevertheless, the enormous nu- 
merical difficulties connected with the calculation of the 
electron-phonon interaction has de-facto prevented the 
systematic application of the HAC theory. Nowadays, 
the advent of more refined numerical techniques, has 
made possible to ground the HAC approach in a fully 
ab initio framework More recently, on this journal, 
Giustino et al. 0, Hcf found a large zero-point renor- 
malization (615 meV) of the band-gap of bulk Diamond. 
The HAC approach is, however, based on a static theor 
of the electron-phonon coupling and, Giustino et al 
in their work, rise a worryingly doubt: We note that the 
good agreement between our calculations and experiment 
may be somewhat fortuitous since the HAC theory does 
not take into account dynamical effects. 

In this work we show, indeed, that when the 
temperature vanishes, the quantistic zero-point mo- 
tion of the atoms induces strong dynamical effects on 
the opto-electronic properties of Diamond and trans- 
polyacetylene, a paradigmatic nano-structure. The sub- 
gap peaks appearing in the experimental absorption spec- 
trum of Diamond are interpreted in terms of polaronic 
states, composed of entangled electron-phonon pairs, 
that cannot be described by the HAC theory. When the 
reduced dimensionality of the system enhances the am- 
plitude of the atomic vibrations, the zero-point motion 
effect even fragments the electrons in a continuum of po- 
laronic states. This is the case of irans-polyacetylene, 
where the single-particle band-structure is replaced by a 
jelly-like electronic distribution and the wave-functions 
of electrons and atoms are stretched along the poly- 
mer axis. This is the break-down of the band-theory. 
By disclosing the physical motivations of the zero-point 
motion effect, we discuss how the present results lead 
to potentially ground-breaking consequences on our un- 
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derstanding of the opto-electronic properties of carbon- 
based nano-structures. 

The zero-point motion effect can be understood by us- 
ing simple arguments. At finite temperature the atoms 
oscillate around their equilibrium positions. These os- 
cillations can be mapped in a system of non-interacting 
harmonic oscillators. As a consequence of the quantistic 
nature of the atoms, when the temperature goes to zero 
the atoms collapse in a ground state with a finite energy. 
This is the zero-point energy that induces quantum fluc- 
tuations of the atoms. The effect of this purely quantis- 
tic motion on the electronic levels can be described using 
MBPT applied to the electron-phonon problem [(I . It 
can be shown that, at zero temperature, the electron- 
phonon self-energy, taken to the lowest non vanishing 
order of perturbation theory, is composed of two contri- 
butions. The first term is the Fan self-energy 
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with | J) and |A) the generic electronic and phononic 
J — f. and = 1 — fj, with fj the electronic 



state. N£ = fj 



occupations. N q is the total number of transferred mo- 
menta used to integrate the self-energy. The key quan- 
tities in Eqll] are the electron-phonon matrix elements 
g$j = (I\Hei-ph\J)g>\\), with H e i- P h the electron 
phonon Hamiltonian. The gjj are calculated ab ini- 



tio using Density- Functional-Perturbation-Theory 11 1 . 



A frequency independent Debye-Waller term, T,j is 
added to Yjf an in order to preserve the translational in- 
variance of the theory. More details about the Debye- 
Waller term can be found, for example, in Ref. (lpj . 

The full frequency-dependent Green's func- 
tion Gi (uS) is readily defined to be Gi(w) = 
(w-£ / -£f a "(w)-£f w/ )~ 1 . The true single- 
particle excitations of the system are obtained as 
poles of Gj. The HAC theory and the more general 
Quasi-Particle approximation (QPA) can be obtained 
making approximations on the frequency dependence of 
T,f an (ui). However we would like to follow a different 
path. 

Physically EqQ] describes the scattering of the bare 
electronic state \I) with the continuum of phonons that 
surrounds the state | J). These scatterings are weighted 
by the coupling terms g t j. This scenario can be described 



using the well-known Fano theory .12] which describes, in 
general, the coupling of a discrete state (|/)) with a set 
of final states surrounded by a continuum of excitations 
( | J } (8> | A ) ) . The scattering with many possible final states 
induces interference effects. When these effects are small 
the term with J = I dominates the sum in Eq[T] and the 
Fano theory predicts the spectral function (SF) Aj (cj) = 
^ |3f [Gi (w)]| to be a Lorentzian. This is the QPA, and 
the center of the Lorentzian is the Quasi-Particle (QP) 
energy, defined to be si + Z / Ef a,l (e / ) + E? w , with 



the rcnormalization fac- 



tor [6j]. If in addition | gfj | <§; ujd, with lod the Debye 
frequency, the Lorentzian width goes to zero, Z = 1 and 
we have the HAC approach. This corresponds to a fre- 
quency independent £ Fcm . When the coupling factors 
g^j are large the interference effects force the system to 
create coherent packets of electron-phonon pairs, 
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We define |P) a polaronic state with energy E P . Indeed, 
it can be demonstrated [l| that, by choosing the coeffi- 
cients Af and Sj A as eigenvectors of an effective Hamil- 
tonian, the Green's function corresponding to Eq[Tl can 
be rewritten as Gi(lu) — J^p | Af | 2 (w - E P + iQ + )~ 1 . It 
also follows that (P | P) = 1. Physically Eq. [2] describes a 
system where the electrons are replaced by a continuum 
of polarons whose purely electronic part is weighted by 
| Af | . The QPA is recovered when the | Af | factors have 
a Lorentzian distribution around the QP energy. How- 
ever, in general, if the \Af\ are small the second term 
in the r.h.s. of Eq.[2]is large making the A](u>) to deviate 
from the simple Lorentzian lincshapc. In this case the 
electron is fragmented in several polaronic states and the 
QPA is expected to fail. 
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FIG. 1. The experimental [3 d?£2 (w) /dw 2 of bulk Diamond 
near the absorption edge (red spheres) is compared with the QPA 
(dashed line) and with the full— dynamical theory (full line). The 
optical gap is defined by transitions between the T' 2 § v (occupied) 
and Tisc (empty) states. The r' 25v state is a genuine QP state. The 
Tl5c SF, instead (shown in the inset), is characterized by strong 
side-peaks that well reproduce the sub-gap peaks appearing in 
the experimental spectrum below 7.1 eV. E^ w represents the QP 
direct gap obtained by neglecting the zero-point motion effect. 

In the inset of FigQ] the SF Aj (w) of bulk Diamond, 
obtained by using the full dynamical dependence of the 
self-energy, is compared with the QPA for the Tisc state, 
the bottom of the conduction bands (l5j. The SF, instead 
of being a simple Lorentzian, shows peaks at 7.25 and 
7.95 eV . These peaks are due to polaronic states |P) 
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each carrying a fraction (given by |A^| ) of the total 
electronic charge. These peaks allow to explain the sub- 
gap states observed experimentally Indeed the ri 5c state 
contributes to the onset of the absorption via transitions 
with the top of the valence bands, represented by the T' 25v 
state. The corresponding contribution to the absorption 
edge will be given by the convolution of the r!,^ and 
]?i5 C SFs. This convolution is used, in Fig.[TJ to calculate 
the second derivative of the dielectric function and to 
compare it with the experimental spectrum [l4| . The 
experimental absorption clearly shows the absorption gap 
at 7.19 eV together with some sub-gap structures at 6.71 
and 6.94 eV. Such structures are completely absent in 
the QPA or using the HAC approach. The dynamical 
theory, instead, leads to an excellent agreement with the 
experimental result. 

This result, while supporting the correctness of a dy- 
namical theory of the zero-point motion effect, ques- 
tions the validity of the QPA and of the static HAC 
approach. The strength of the dynamical zero-point mo- 
tion effect in Diamond is unexpected, and opens the path 
to potentially stronger effects in carbon-based nano- 
structures. Indeed, the strength of the electron-phonon 
self-energy is linked to the amplitude of the atomic vi- 
brations. In a nano-structure the atoms are more free 
to oscillate, thanks to the reduced symmetry of the sys- 
tem. This can be easily verified by considering a con- 
jugated polymer, as irans-polyacetylene. This is a ID 
chain made up of repeated structural units linked by al- 
ternated single and double bonds between the carbon 
atoms 16]. The unit cell contains 2 carbon and 2 hy- 



drogen atoms, which lay on the same plane, as shown 
in the upper frames of FigJ5] We can associate an av- 
erage quantistic size to the atoms by using the standard 
deviation Oi of the atom i in its ground-state wavefunc- 



tion [7]. This is Si « ^Ea^a 1 6 ( A ) > with 6 ( A ) 
the polarization vector of the phonon mode A. In dia- 
mond ac ~ 0.1 a. u., independently on the direction. In 
irans-polyacetylene the smaller distance between carbon 
atoms, slightly reduces the standard deviation along the 
x direction where ac ~ 0.08 a. u. while in the y direction 
ere ~ 0.16 a.u.. Hydrogen is twelve times lighter than 
carbon. As a consequence its standard deviation in the 
polymer plane is an ~ 0.2 a.u.. 

To clearly visualize the dramatic effect of the 
zero-point motion effect on the electronic structure of 
irans-polyacetylene we define a global SF A(k, Li) = 
En w 1^ [Gnk where we have explicitly expressed 

the index / in terms of the band (n) and k-point (fc) in- 
dexes, k is taken in the V ~ X direction [15[. Physically 
AZ = A(k,uj) Aw gives the fraction of electronic charge 
carried by the state with k-point k in the small energy 
range Aw. From the definition of A (k,u>) it follows that 
AZ cx \Af | . The A function is a very peaked function of 
k and w both in the HAC approach and in the QPA with 



the charge confined in very sharp single-particle states. 
The plot of the A function, shown in Figj2]in the energy 
range of the last three occupied bands, gives instead a 
completely different picture. The intensity reported in 
the color scale of the FigJ5] refers to the dimensionless 
quantity AZ calculated with Auj = 50 meV. As a refer- 
ence we show the range of values of AZ corresponding to 
the QPA. Note that the static HAC theory corresponds 
to AZ = 1. 




FIG. 2. Two-dimensional plot of the SF A(k,ui) in the energy 
region of the last three occupied bands. The range of values of A 
are given in terms of the dimensionless quantity AZ (see text) that 
measures the elemental charge carried by the SF. Except the 7r/-7r* 
states, perpendicular to the polymer plane, the electronic levels ac- 
quire a large energy indetermination that is particularly dramatic in 
the region enclosed by the two dashed line. This destructive effect 
is, instead, a cooperative interaction between atoms and electrons 
that is reflected in the stretching of the polaronic states in the di- 
rection of the polymer axis. This is evident from the comparison 
of the wavefunctions of a reference bare state (left frame, point a) 
with two polaronic states (central and right frame, points b and 
c). The standard deviation of the atoms in their unperturbed state 
and in the polaronic states is represented by the size of the spheres 
in the three upper frames. 

With the exception of the last occupied band, the 
charge carried by the SF is lower then 0.4, well below 
the range of values where the QPA or the HAC theory 
is accurate. More importantly very small values of AZ 
mean small \Af\ and, consequently, large contributions 
of the electron-phonon pairs in the polaronic state (sec- 
ond term in the r.h.s of EqJ5]). This effect is even more 
dramatic in the central energy region, where the bands 
disappear being replaced by an almost uniform, jelly-like 
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charge distribution that, in Fig{2j is represented by the 
region enclosed by the two dashed lines. Physically the 
system is forced to split the electronic levels in several 
polaronic states. Single-particle states are not allowed 
anymore. This is an unambiguous signature of the break- 
down of the band theory 

With respect to more deeper states the last occupied 
bands, that has a tt/tt* character, corresponds to states 
distributed perpendicular to the polymer axis. As a con- 
sequence they feel less the effect of the on-plane vibra- 
tions and almost all electronic charge is confined in sharp 
quasi-particle-like peaks. 

A more careful analysis reveals that the destructive 
effect of the zero-point motion effect on the electronic 
band-structure is, instead, a signature of a cooperative 
electron-atom dynamics. Indeed, EqJ5] allows to calcu- 
late, ab initio, several properties of the polaronic states. 
By projecting \P) on the real-space we can define the 
polaronic wave-function, (r\P), and also calculate the 
standard deviation af of the atom i that participates 
in the polaronic state. This can be used to explain the 
physical properties of the wide region enclosed by the two 
das hed linoa in Fig.[2l In the upper frame (a) of Fig. [2] 



particles. This inevitably leads to the failure of the band- 
theory. The present results highlight the limitations of 
an approach to carbon-based nano-structures that ne- 
glects the effect of the electron-phonon coupling. The 
deep change of the electronic and atomic wave-functions 
induced by the formation of the polaronic states can 
have important implications on related properties. This 
inevitably, imposes a critical revision of the results ob- 
tained using purely electronic theories with atoms frozen 
in their equilibrium positions. 
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tional collaboration project. 



we conoidor a bare electronic otato that io localized on 
the C — H bond. The zero-point motion effect splits the 
charge of this state in several polaronic states. Two of 
the most intense arc the states (b) and (c) whose wave- 



fun qtions are shown in the upper frames of Fig. [2] By 



moving towards the center of the region the electronic 



chan ge is gradually moved from the C — II bond towards 



the polymer axis. — In the state (c) thia charge transfer 



is mpre pronounced and the wave-function is completely 
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delocahzed along the polymer axis. Even more stunning 
is the effect of the zero-point motion on the atomic stan- 
dard deviation. In both states (b) and (c) Sf is decreased 
by 50% in the y direction and increased by 150% in the 
direction of the axis. The present results show that the 
zero-point motion effect works to enhance the dereal- 



ization of the charge carriers by stretching both electrons 



and atoms along the polymer axis. This cooperative ef- no 



feet agrees with the general interpretation of the high 



electronic mobility in polymers as dictated by the effi- 
cient hopping of the charge carriers [II El- 



In conclusion, we have shown that, when the temper — 
atur c vanishes, the quantistic zero point motion of the 

atoms induces large corrections to the opto-electronic 
properties of Diamond and frans-polyacctylene. The 
electronic bands are replaced by a continuum of po- 
laronic states composed by entangled electron-phonon 
pairs. This successfully explains the sub-gap states ob- 
served experimentally in the absorption spectrum of Dia- 
mond, and the rich structures that appear in the spectral 
functions of irans-polyacetylene. The cooperative dy- 
namics between electrons and atoms that leads to the for- 
mation of the polaronic states rules out any description 
in terms of bare atoms, bare electronic states or quasi- 
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